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Abstract 

A three-dimensional chemical transport model (CTM) 
(Kaminski, 1992) has been used to study the evolution of 
the Arctic ozone during the winter of 1992. The continu- 
ity equation has been solved using a spectral method with 
Rhomboidal 15 (R15) truncation and leap-frog time step- 
ping. Six-hourly meteorological fields from the Canadian 
Meteorological Centre global objective analysis routines 
run at T79 were degraded to the model resolution. In ad- 
dition, they were interpolated to the model time grid and 
were used to drive the model from the surface to 10 mb. 
In the model, processing of Cl r occurred over Arctic lati- 
tudes but some of the initial products were still present by 
mid-January. Also, the large amounts of CIO formed in 
the model in early January were converted to CINO3. The 
results suggest that the model resolution may be insuffi- 
cient to resolve the details of the Arctic transport during 
this time period. In particular, the wind field does not 
move the C10 r ‘cloud' to the south over Europe as seen in 
the MLS measurements. 


Introduction 

Although the dynamical situation during the Arctic 
winter is less stable than that in Antarctica and the tem- 
peratures are not as cold, Polar Stratospheric Clouds (PSC, 
type I and type II PSCs do form (eg. GRL, 1990). Also, air 
remains in the vortex sufficiently long that there may be 
extensive processing (eg., McConnell et ah, 199f) so that 
Cl r may be converted to active forms of chlorine within 
extensive regions of the polar vortex. This has clearly 
been observed during the Arctic Airborne Stratospheric 
Expedition (AASE) during the winter of 1989/90 where 
enhanced levels of HrO and (TO were observed along with 
some denitrification (GRL, 1990). Further, the microwave 
limb sounder (MLS) instrument on board the EARS satel- 
lite has observed extensive ‘clouds’ of (JO over northern 
latitudes and as far south as the Mediterranean (Waters 
et ah, 1992). 


In this note we present the results of an attempt to sim- 
ulate the dynamics and chemistry of Arctic stratosphere, 
and in particular the formation of t he CIO ‘cloud’ observed 
over Northern Europe during the 1991/1992 winter (eg. 
Waters et ah, 1992). The York Chemical Transport Model 
(CTM) (Kaminski, 1992) was run in an off-line mode from 
December 15th to January 22nd. The CTM was driven 
by the objectively analyzed wind and temperature data 
generated by the Canadian Meteorological Centre (CMC) 
global spectral model (Ritchie, 1991) run in the Analysis 
Cycle (zero hour forecast at G0Z, 06Z, 12Z, and 18Z was 
used). The CMC’s global objective analysis (OA) stops 
at 10 mb. The upper level temperature fields are inverted 
from satellite observations (TOVS), and the wind fields 
are derived from temperature. 

The chemistry in the model consists of 49 reactions 
(which includes gas phase, heterogeneous and photochem- 
ical reactions) and 21 species. The following heteroge- 
neous reactions were included 

C10N0 2 (g) + HCl(s) - Cl 2 + IIN0 3 (s) (1) 
N 2 0 5 (g) + HCl(s) CIN0 2 (g) + HN0 3 (s) (2) 

C10N0 2 (g) + H 2 0(s) - HOCl(g) + HN0 3 (s) (3) 

N 2 O s (g) + II 2 0(s) -+ 2HN0 3 (s) (4) 

and were activated by temperature switches. For T < 195 
K the reactions (1) and (2) were activated, while for 
T < 190 K the reactions (3) and (4) were added. Al- 
though this does not simulate the details of the effects 
of processing by type I and II PSCs it should produce 
the correct general conversion from inactive forms of Cl x 
to active forms, CIO*. We define Cl r = (TO x + 1 1 Cl T 
CINO3 and C10 x = Cl + 2C1 2 + CIO + IIOC1 + C1N0 2 
T 2C1 2 0 2 . 

Seven tracers are used: O x , Cl x , N0 2 , (TI.|, N 2 0, H 2 0 
and passive- 0 3 . The passive tracer lias the mixing ratios 
of the initial 0 3 field and provides a. measure of the ef- 
fects of transport without chemistry and is thus particu- 
larly useful for comparing with the chemically active 0 3 
field. In addition, it provides a useful check for the mass 
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conservation properties of the code. For this experiment 
a half-hour dynamical time step (required by the CFL 
condition) and a 1.5 hour chemical time step are used to 
calculate the species concentrations. The chemical time 
step of 1.5 hours was selected in order to adequately sim- 
ulate the diurnal behaviour of the Arctic chemistry. The 
chemistry solver is stable and previous experiments have 
indicated that much longer time steps could be used. The 
J values are evaluated by means of a look-up table at 
each chemical time step. A two-dimensional look-up ta- 
ble is generated using combinations of the logarithm of O 2 
and O 3 optical thickness, and temperature, depending on 
the J value. Currently we allow the J value look-up ta- 
bles to have a maximum error of about 5 % as compared 
with those generated by a detailed wavelength integration 
method (cf. Henderson et ah, 1987). 

The CTM employs a spectral solution to the conti- 
nuity equation (eg. Williamson et al., 1987) and was 
run with Rhomboidal 15 (R15) truncation, correspond- 
ing to 48 longitudes by 40 latitudes in grid space. A 
leapfrog scheme in spectral space, and a finite difference 
scheme in the vertical are used to integrate the continuity 
equation. The initial conditions, representative for mid- 
December, were generated by a two-dimensional zonal 
average model of the middle atmosphere that was sup- 
plied by Dr G. Brasseur of NCAR (Brasseur et al. 1990). 
The model was run from 15th December, 1991 to 22nd 
January, 1992. Since this is a short note and the results 
are extensive we have chosen to concentrate on the out- 
put for a single day and a single pressure level, viz. the 
30 mb surface on the 16th January, 1992. 


Analysis and Conclusions 

It is quite clear from a preliminary comparison of the 
model O 3 data with the TOMS 0 3 that the dynamical 
field rapidly imposes its main features in a relatively short 
period. A comparison between the model active and pas- 
sive ozone and the TOMS ozone field is shown in Fig- 
ure 1 for the 16th January, 1992. It can be seen that 
many of the features exhibited by the TOMS ozone are 
reproduced by the model. Noteworthy are the highs over 
eastern North America and eastern Asia and the low over 
the eastern Atlantic. There are several locations where 
the detail is not good. However, considering the limited 
resolution of the model and the fact that it presents a 
snapshot whereas the TOMS data represent the ozone 
field gathered over a day we consider the agreement quite 
reasonable. Since the effects of chemistry in the polar re- 
gions are marginal over this short time period what we 
are seeing in the model is largely the effects of transport. 
This is also seen in a comparison of the active and passive 
ozone fields in Figure 1 . At low latitudes the difference 
between the passive and active ozone is due to the ef- 
fect of substantial ozone generation, as expected, due to 


photolysis of 0 2 . The effect of ozone depletion at boreal 
latitudes is not immediately evident and still requires fur- 
ther quantification. We note that the model ozone highs 
are within 5% of the TOMS data. 

At the beginning of the simulation the temperature 
is sufficiently cold to form PSC type I and the hetero- 
geneous reactions listed above are activated. During the 
second week of January, 1992 the boreal temperatures are 
warmer than 195 K and there is no further formation of 
active chlorine. Figure 2 shows time series of the percent 
volume of air within each layer that has (a) temperatures 
less than 190 K, (b) temperatures less than 195 K and (c) 
temperatures less than 200 K for latitudes greater than 
50° N. It shows quite dramatically the evolution of the 
processing in the model. The 200 K temperature was in- 
cluded because the CMC temperatures in this region of 
the stratosphere for this period of the forecast were prob- 
ably up to 5 degrees too warm (Dr H. Ritchie, private 
communication, 1992; this problem with the stratospheric 
temperatures has now been corrected). 

Although the processing has finished by the 16th of 
January it takes several weeks for the processed air to 
be reconverted back to HC1 and CINO 3 (eg. Henderson 
et al., 1990). This is indicated in Figure 3 which shows 
the fields for HC1, C10 r , and CINO 3 on the 30 mb level 
for OOZ on the 16th January, 1992. Figure 3a shows that 
a large region of HC1 has been completely processed to 
active chlorine and has not been reformed. At the time of 
the snapshot the hole in the HCl field is over Greenland 
and northern Europe. Other regions with low values are 
due to transport effects which are immediately apparent 
in the Cl x field (not shown). Furthermore, we note that 
there as almost no PCS type II activated in the model. 

The active chlorine field, C10 x , is shown in Figure 3b. 
It is similar in shape to the processed HCl hole, but some- 
what smaller and this is significant as illustrated later. 
The peak C10 x mixing ratio is 3.6 ppb. The CIO* >• 
composed mainly of C1N0 2 , HOC1 and Cl 2 with peak 
mixing ratios ~ 2.2 ppb (over Greenland), ~ 0.45 ppb 
over Scandinavia and ~ 0.6 ppb over Greenland, respec- 
tively. C1 2 0 2 and CIO are relatively small components 
of the C10 x field. The fact that C1N0 2 , HOC1 and Cl 2 
predominate, indicates that the air in this region has not 
seen the sunlight otherwise it would have been processed 
into CIO and C1 2 0 2 . Given that the MLS measurements 
of Waters et al. (1992) show several ppb of CIO, our re- 
sults indicate that there may be serious problems with 
the dynamical field run at R15. 

The missing Cl x between the HCl hole and the CIO* 
‘cloud’ is found as CINO3 and is shown in Figure 3 c. It 
appears as a collar of C 1 N 0 3 around the CIO* ‘cloud’ with 
the highest values ~ 3 ppb occurring over northern Russia 
and with mixing ratios ~ 1 ppb around the collar. What 
appears to have occurred is that HCl has been processed 
and the products photolyzed to form CIO. This has then 
reacted with the N 0 2 present in the model to form CINO3. 
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Figure 1 . a) Active, b) passive and c) TOMS column ozone for the 16th 
January, 1992 
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Figure 2. Time series showing the percent volume of air that has a) temper- 
atures less than 190 K, b) temperatures less than 195 K and c) temperatures 
less than 200 K for latitudes greater than 50° N (time sries are shown for 4 
model levels). 
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Figure 3. a) The HC1 mixing ratios on the 30 mb surface for 00Z on the 16th 
January, 1992 using the CTM showing the region of reduced mixing ratios 
due to processing simulated by the model; b) ClO x c) CINO 3 . 






We note that we have not included any processing of N 2 0 5 
in the model to HNO 3 by stratospheric aerosols which 
would occur at all temperatures. This would remove any 
N 2 O 5 formed in the model and so the N0 X . We find that 
N 2 O 5 mixing ratios are ~ 1-3 ppb in the polar region 
(> 50° N) and NO x values are < 2.5 ppb. However, this 
would be adequate to convert any free CIO to CINO 3 
on a time scale of a day. The edge of the CINO 3 collar 
reaches as far south as northern France. This is about 
1 model grid point further north than measured by the 
MLS instrument for earlier dates. Again, we suspect that 
this is due to the limitations of the R15 simulation and 
possibly poor representation of the N 2 O 5 as mentioned 
above. 
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